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A single-column model run under the weak temperature gradient 

' ' 9 approximation, a parameterization of large-scale dynamics appro- 

\ 10 priate for the tropical atmosphere, is shown to have multiple stable 

11 equilibria. Under conditions permitting persistent deep convection, 

O I 12 the model has a statistically steady state in which such convection 

3 occurs, as well as an extremely dry state in which convection does 

4 not occur. Which state is reached depends on the initial moisture 

5 profile. 
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1. Introduction eo 

81 

In the extratropics, balanced dynamics associated withi barOg^ 
clinic eddies can force precipitation strongly by inducing ascer^^ 
and adiabatic cooling. In the tropics, the reasons for the onset o^^ 
precipitation at a given time and place are often much more subtle^^ 
and large-scale ascent more a response to deep convection than a 
cause. On the other hand, in the tropics the relationship of time- 
mean precipitation to boundary conditions is generally stronger 
than it is in the extratropics. In simple models of tropical climate 
dynamics, it is typical to treat the time-mean tropical precipitation'' 
as deterministically related to the boundary conditions. While this" 
may for many purposes be adequate, it is clear that the time mean' 
state, particularly in climatologically rainy regions, is an averagl^ 
over periods of strong precipitation and periods of negligible pre-^ 
cipitation. As precipitation cannot be negative, this implies an in-'' 
herent nonlinearity that is overlooked in deterministic models o¥ 
the time mean state. It also suggests that transient disturbances, o?^ 
other factors whose relationship to the boundary conditions is aft'' 
best hard to discern, may play a role in maintaining the time-meaff 
state. ^ 

In this study we present a very simple expression of these prop- 
erties of the tropical atmosphere, namely its nonlinearity and the 
potential complexities hidden in its time-mean response to bound- 
ary conditions. We show that a single-column model with boundary 
conditions permitting strong deep convection can have two stable 
steady (or statistically steady) states: one in which persistent deeg^ 
convection occurs, and one in which it does not. Which solution 
is reached in numerical integrations depends on the initial moisture^ 
field. 

104 

105 

2. Model and Experiment Design loe 

107 

We use the single column version of the GEOS-5 genos 
eral circulation model (GCM). A description of the GEOS09 
5 system and its physical parameterizations may be found ^ 
http://gmao.gsfc.nasa.gov/systems/geos5/, and a brief outline ift 
presented in the auxiliary materials. 112 

The model is nm over an ocean surface with fixed SST. The iijifj 
solation and solar zenith angle are held constant at values of 4Qf)^ 
Wm~^ and zero respectively. The surface wind speed, used in tfj^ 
bulk formulae for the surface fluxes, is set to a constant 7ms~^. 

The experiment design is essentially the same as that used i,^ 
5obel and Bretherton [2000; hereafter SBOO], The model is firfi^tj 
run to a state of radiative-convective equilibrium (RCE) over a^r^^ 
SST of 301 K. In this calculation the large-scale vertical velocit^^ 
is set to zero. The temperature profile from the RCE calculation is 

then used as an input to a set of calculations in which the model is 

122 

modified to implement the weak temperature gradient (WTO) ajii- 
proximation. In the WTG calculations, the temperature profile is 
held fixed in time in the free troposphere, defined somewhat arbi-'' 
trarily as those levels with pressures less than pinv= 850 hPa. A¥ 
those levels, the large-scale vertical velocity is diagnosed as thJif 
which causes the vertical advection of potential temperature to prh^- 
cisely balance the diabatic heating computed by the model physic^,^ 
consistent with the requirement of zero temperature tendency. Ko)?- 
izontal advection of temperature is also assumed negligible. In tHS" 
nominal boundary layer, defined as levels with pressures greatii' 
than pinv, the temperature is determined prognostically, with thS^ 
vertical velocity computed by linear interpolation in pressure 
tween the diagnosed value at pinv and an assumed value of zero SI* 
the surface. At all levels, large-scale vertical advection of humidity 
and condensate — both of which are integrated prognostically ass 
usual — are computed using the vertical gradients derived from the? 
internally predicted profiles of those variables and the large-scales 
vertical velocity diagnosed as described above. 139 
In some experiments, the initial moisture profile is taken froiao 
the steady-state RCE solution. In others, the initial moisture profilei 
is obtained by setting the free-tropospheric humidity to zero. 142 



In some experiments, horizontal advection of moisture is ne- 
glected, as in SBOO. In others, it is parameterized by a relaxation 
back to a fixed profile - the steady-state RCE solution - with a fixed 
time scale. The time scale can be thought of as an advective one 
given by a length scale associated with horizontal moisture gra- 
dients divided by a velocity, where the velocity is thought of as 
a rotational one which is independent of height (so that the time 
scale itself is independent of height), and also independent of the 
magnitude of the divergent circulation implied by the vertical ve- 
locity in the column. Another approach is to consider advection by 
that divergent circulation itself, in which case the relaxation rate 
is simply dui/dp, where cj is the pressure vertical velocity and p 
the pressure [Raymond and Zeng, 2005]. We have performed cal- 
culations with this method also; the results are not quaUtatively 
different from those using the fixed advective time scale. 

All simulations are run for one year, by which time the model 
has reached a solution that is statistically steady, and in many 
cases close to truly steady; oscillations sometimes occur, as de- 
scribed briefiy below. Results presented are averages over the last 
2 months. 

3. Results 

3.1. GEOS 5 SCM 

Fig. 1 shows the temperature and relative humidity profiles ob- 
tained from the RCE calculation. The temperature profile is close to 
moist adiabatic as expected. Fig. 2 shows time-mean precipitation 
as a function of SST for the WTG experiments. The figure is simi- 
lar to fig. 4 of SBOO; each point represents a different experiment, 
each of which uses the same free-tropospheric temperature profile, 
obtained from the RCE calculation. The difference here [besides 
that we use a different model; SBOO used the model of Renno et al. 
[1994]] is that we show several sets of curves, obtained using differ- 
ent initial moisture profiles, and different time scales for horizontal 
moisture advection (SBOO did not include this process in any of 
their calculations). We first describe the solid curves, which were 
obtained using no horizontal moisture advection. 

When the initial moisture profile is taken from the RCE, a rainy 
state is obtained for sufficiently large SST. The shape of the curve 
of precipitation vs. SST for this set of solutions (the upper blue 
curve) is qualitatively similar to that obtained by SBOO. The lower 
blue curve curve shows solutions obtained using an initial moisture 
profile is set to zero in the free troposphere. 

When the initial profile is sufficiently dry, apparently a dry state 
can be maintained although the SST is high and the convective 
available potential energy (CAPE) for an imdilute parcel ascent 
would be substantial. For a dry profile, the convective parameter- 
ization is unable to generate significant convective heating or pre- 
cipitation, presumably due to inhibition by entrainment of the dry 
air. In the absence of convective heating, radiative cooling forces 
descent, which maintains the dry state in the free troposphere. On 
the other hand if the initial conditions are sufficiently moist, for suf- 
ficiently large SST the deep convective scheme is able to become 
and remain active. If the ensuing heating exceeds the radiative cool- 
ing, it induces large-scale ascent, which moistens the atmosphere 
further by large-scale advection, leading to the maintenance of the 
convective state. For SST below a critical value, here around 300.5 
K, there is insufficient CAPE (or strictly, work function; Arakawa 
and Schubert 1974) for deep convection to occur even if the initial 
profile is moist. Dry air can suppress deep convection in the pres- 
ence of large CAPE, but moist air cannot cause it to occur if CAPE 
is absent. Thus the multiple equilibria occur only for SST above 
the critical value. 

Fig. 3 shows vertical profiles of relative himridity as functions of 
SST in the dry state and the moist state, without horizontal moisture 
advection. We see that the relative humidity in the boundary layer 
is similar in both sets of solutions, but that the free-tropospheric 
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143 humidity is essentially zero in tlie dry case; in tlie absence of lro»3 

144 izontal moisture advection, tliere is no source of moisture to bal»4 

145 ance drying due to descent. This is somewhat artificial, as in thes 

146 real atmosphere even the driest desert is connected to regions of fiss 

147 nite humidity by horizontal advection. The magenta and red curves 

148 in fig. 2 show results using horizontal advective time scales of 3 

149 and 6 days, respectively (the no-advection case would correspond 

150 to an infinite advective time scale); again the relaxation is towards 

151 a target moisture profile equal to that in the RCE. We see that for 

152 the 3-day time scale, multiple equilibria exist only over a narrow? 

153 range of SSTs. The horizontal advective moistening in this casjea 

154 prevents the atmosphere from becoming sufficiently dry to inhibit 

155 deep convection once SST exceeds the RCE value of 301 K by a 

156 small increment. On the other hand, for a mixing time scale of 6 

157 days, the multiple equilibria persist up to large SST. In this case, 

158 the horizontal advective moistening is too weak to overwhelm the 

159 vertical advective drying, and the atmosphere can stay dry enough 

160 to inhibit deep convection even for large SST if it is initialized suf- 

161 ficiently dry. These results show that the existence of multiple equij^j 

162 libria is not purely an artifact of the neglect of horizontal moistu^g^ 

163 advection, although sufficiently strong horizontal moisture advegjj 

164 tion can eliminate the dry equilibrium (or render it unstable so that 

165 it cannot be achieved numerically). 

166 In the cases with horizontal moisture advection, the precipita- 

167 tion is a non-monotonic function of SST, with local maxima near 

168 301K, the value at which the RCE was computed. The solutions 

169 in this neighborhood are time-dependent, either periodic or chaotic 

170 but with a dominant spectral peak. The periods range from 1-20 h, 

171 and maximum excursions of 0.5-10 mrndr^ , with the largest os- 

172 dilations tending to occur for the higher SST values although the 

173 dependence of amplitude on SST is not monotonic. We do not un- 

174 derstand these oscillations in detail and do not address them further 

175 here; our interest is in the presence or absence of multiple equiU^gj 

176 ria. 



3.2. Interpretation 
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177 We consider under what conditions the multiple equilibria caff 

178 exist. We can derive some general constraints that any model musP 

179 obey if it is to obtain a dry solution under boundary conditions an3^ 

180 forcings that also allow a rainy solution. 

Starting from the primitive equations, the steady temperature 
and moisture equations in pressure coordinates are 
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u • VT ojdps = Q^ + R- dpoj's', 
u • Vg -I- udpQ = Qq — dpui'q', 



214 
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181 where u is horizontal velocity, T temperature in energy units (i.es.t? 

182 multiplied by the heat capacity of air at constant pressure Cp), s djys 

183 static energy, q specific humidity in energy units (i.e., multiplisio 

184 by the latent heat of vaporization L„), R radiative heating, Qc cofl^i 

185 vective heating, Qq convective moistening, uj large-scale vertic# 

186 velocity, and V the horizontal gradient on pressure surfaces, aj'??^ 

187 and uj'q' are the turbulent fluxes, limited to the atmospheric bound^'' 

188 ary layer (ABL), as parameterized by the model's boundary laySi* 

189 scheme; in the free troposphere these fluxes are incorporated in6S'' 

190 QcandQq. 

The WTG approximation requires neglect of horizontal temp#^ 
ature advection in the free troposphere; for the sake of argume^ 
let us also neglect horizontal moisture advection and horizontM" 
temperature advection in the ABL. If there is no deep convection^' 
Qc = Qq = 0, the equations for the free troposphere are: 



ujdpS = R, 
ojdpq = 0. 



232 
233 
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191 Therefore, unless the radiative cooling R is zero, the vertical vese 

192 locity u) is non-zero and dpq = 0. Integrating from the top of thss 



atmosphere, this in turn yields q = : the free troposphere has 
to be dry for the non-convective equilibrium to exist. If horizontal 
moisture advection is included, the free-tropospheric humidity can 
be non-zero in this dry state. 

Equations (3) and (4) can be integrated over the ABL: 



{ujdps) = {R) + H, 

{u)dpq) = E, 



(5) 
(6) 



where H and E are the surface sensible and latent heat fluxes, and 
indicates the integral from the surface to the top of the inversion 

Pinv. 

In our formulation, u varies linearly from pinv to the surface. 
Equations (5) and (6) can therefore be rewritten: 



WinvISM - Sinv) = + H, 

l^inv qM = E, 



(7) 
(8) 



where the subscript M indicates the vertically averaged value over 
the ABL and is the dry static energy at the top of the inversion. 
We used the fact that the humidity at that level is zero. 

Equation (8) gives us one constraint that must be satisfied in or- 
der for the dry solution to occur: the ventilation of the ABL by the 
subsidence has to compensate the evaporation. On the other hand. 
Equation (7) shows that the warming by surface sensible heat flux 
and subsidence has to be compensated by the radiative cooling. The 
two equations can be rewritten as a single condition for the exis- 
tence of the non-convective equilibrium: 



{R)+H 



qM 



E. 



(9) 



As the SST increases, sm increases, while is fixed in WTG. 
The ABL stability — sm thus decreases. The radiative cooling 
— {R) increases too, and the surface heat flux H is not very sensi- 
tive to SST changes. The subsidence Winv therefore must increase 
with increasing SST. The ABL specific humidity qM also increases 
with SST. So the effect of ventilation increases with SST. 

To maintain the dry solution over an ocean surface, it is impor- 
tant that the ABL air be able to stay relatively moist as the SST 
increases so that E is limited and the ventilation a;inv?M is effi- 
cient. Otherwise E will become very large due to the large air-sea 
humidity contrast (taking surface wind fixed), and in general there 
is no mechanism for R nor (s[^^ — sm)~^ to become large at the 
same time to compensate. In nature (and in many models), the 
ABL generally does stay moist under a dry free troposphere, e.g., 
in subtropical trade wind regions. In some models, though, this 
may not occur; obvious examples are idealized models in which the 
shape of the entire vertical structure of the himiidity field is fixed, 
so that the ABL moisture is proportional to the free-tropospheric 
moisture. This happens, for example, in the first quasi-equilibrium 
tropical circulation model (QTCM; Neelin and Zeng [2000] ; Zeng 
et al. [2000]). This can be remedied by allowing a separate degree 
of freedom for boundary layer moisture, as in some otherwise sim- 
ilar models [Wang and Li, 1993; Meggers et al., 2006; Sobel and 
Neelin, 2006; Khouider and Majda, 2006]. 

The second constraint that must hold in order for the dry solu- 
tion to occur is that the convective heating and moistening remain 
zero, or at least small. For a dry free troposphere and a cold sea sur- 
face, so that CAPE is negative, this will occur for any reasonable 
convective parameterization. For a dry free troposphere and a warm 
sea surface, so that non-entraining CAPE is significantly positive, 
the results may be model-dependent. Parameterizations which are 
insufficiently sensitive to dry free tropospheric air may be able to 
generate some heating. Once heating occurs, ii Qc + R > there 
will be ascent, which will moisten the troposphere, leading to the 
establishment of the rainy solution. Alternatively, if the scheme is 
not able to generate deep convection and associated heating, but 
can produce enough moistening above the ABL to eventually allow 
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240 deep convection to occur, this will also lead to the rainy solutioais 

241 It is reasonable to expect that under the same boundary conditions^ 

242 some models will have multiple equilibria and others will not. 

243 We have performed preliminary experiments with a small num- 

244 ber of other SCMs in addition to the GEOS5. In these experiments, 

245 we have been able to produce multiple equilibria in some but not 

246 all of the models. However, our search over the space of initial 

247 and boundary conditions in these other models has been cursory," 

248 so it is premature to draw any conclusions from these results. Wfe^ 

249 will report on them in more detail when we have performed a moii? 

250 thorough study. 

314 
315 

4. Conclusions 316 

317 

251 We have shown that a single-column inodel using essentially thes 

252 same physics and numerics as a state-of-the-art GCM has multiples 

253 equilibria, when run with fixed free tropospheric temperature arsdo 

254 diagnostic large-scale vertical velocity according to the weak tem^i 

255 perature gradient approximation. When the boundary conditiofl82 

256 are such as to allow a rainy equilibrium state, a second equilibrium' 

257 with no precipitation and a very dry free troposphere also exists,"* 

258 and is reached by initializing the model with a very dry soundings 

259 This dry equilibrium state persists in the presence of paramete!^^ 

260 ized horizontal moisture advection, represented as a relaxation &f 

261 the specific humidity back to a relatively moist reference profile, If 

262 long as the relaxation time scale is not too short. The critical valiie" 

263 of this time scale for existence of multiple equlibria is about four" 

264 days in this model. When the SST is sufficiently low, compared to_ 

265 that at which the free-tropospheric temperature profile would be 

266 radiative-convective equilibrium (RCE), only the dry equilibriugj, 

267 exists. 335 

268 The existence of the dry equilibrium under SST greater th;^ 

269 the RCE value requires that the radiatively-driven large-scale d§jj 

270 scent be able to export sufficient moisture to balance the surfaQ§8 

271 evaporation, and that the convective parameterization (or explicife 

272 convection, in a cloud-resolving model) be sufficiently sensitive to 

273 free-tropospheric moisture that the dry troposphere inhibits deefn 

274 convection from occurring. Whether these requirements are me* 

275 may be model-dependent. 313 

276 The existence of these multiple equilibria is a consequence 

277 the interaction between deep convection and large-scale dynamics*,^ 

278 with the latter parameterized here through the weak temperatu?6^ 

279 gradient approximation. The multiple equilibria presented here aiff 

280 in this respect fundamentally different than those found by Renii^ 

281 [1997], whose model did not include a representation of large- scafS" 

282 dynamics. It seems possible, but is not obvious, that analogoifl" 

283 multiple equilibria would exist for other parameterizations of lar^e^ 

284 scale dynamics, similar in spirit to the WTG approach used here b^|j 

285 differing in detail (e.g., Bergman and Sardeshmukh [2004]; Mapes^ 

286 [2004]; Kuang [2007]). 3^5 

287 The existence of these multiple equilibria strikes us as a direg^ 

288 and simple expression of the tropical atmosphere's inherent nonligjy 

289 earity. It illustrates the complexity that is hidden behind the averag^a 

290 ing when the time-mean precipitation, for example, is consider^ 

291 as a function of boundary conditions. 

292 We cannot be certain that the multiple equilibria would exist 

293 for a hypothetical model with "correct" physics. Recent studies 

294 suggest, however, that the tendency of many current convective pa- 

295 rameterizations is to have too little, rather than too much, sensitivity 

296 to free tropospheric moisture, due to insufficient entrainment (e.g., 

297 Derbyshire et al. [2004]; Kuang and Bretherton [2006]; Biasutti 

298 et al. [2006]). As we expect that more sensitivity will make the 

299 multiple equilibria more likely to exist, this suggests that they are 

300 not purely an artifact of a bias in the GE0S5 physics. In addition, 

301 in recent simulations of radiative-convective equilibrium on a large 

302 domain over uniform SST [Bretherton et al. , 2005] deep convection 

303 occurs only in a single, small region, while the rest of the domain 

304 becomes extremely dry. This behavior is analogous to what we find 

305 here, with our single column representing either the dry or rainy re- 

306 gion separately. This suggests that a cloud-resolving model run in 

307 WTG mode might well exhibit multiple equilibria. 
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Figure 1. Temperature (a) and relative humidity (b) profiles in 
the radiative-convective equiUbrium state over an SST of 301K. 
Temperature is in °K, relative humidity in %. 
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Figure 2. Precipitation (mmd~^) as a function of SST in tlie 
convective and non-convective equilibria, witli no horizontal 
moisture advection (blue, solid) and with horizontal moisture 
advection parameterized as a relaxation back to the RCE profile 
with a time scale of 6 days (magenta, dashed) and 3 days (red, 
dot-dashed). Each circle represents a separate model integra- 
tion. 
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Figure 3. Relative humidity (%) as a function of SST and pres- 
sure in tlie non-convective (a) and convective (b) equilibria, with 
no horizontal moisture advection. 



